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QUATERNION DEVELOPMENTS WITH APPLICATIONS* 

BY 

JAMES BYRNIE SHAW 

Introduction. 

This paper is an extension of the quaternion algebra along lines analogous 
to those of a preceding paper, f Certain functions of one or more quaternions 
are studied in detail, all of them independent of the orientation of the unit 
axes, 1, i,j, k. The notation is functional and seems preferable to the use 
of brackets such as those used by Hamilton, J who, indeed, called his 
notation temporary. We seem to have here a natural method for treating 
several functions that occur frequently in analysis. For example, if we 
represent quaternion numbers by vectors in four-dimensional space, as 
was done by Stringham§ and Hathaway, || these functions are expressions 
for geometrical relations which are independent of the particular system of 
rectangular axes used. These expressions are useful in electrodynamics and 
in relativity problems. They include the Hamiltonian [ ab ] , ( abc ) , [ abc ] , 
(abed) , and the equivalent forms used by JoLY,1f and Schrutka von Recht- 
enstamm.** 

The development of these formula? enable us also to formulate certain dif- 
ferential operators of frequent occurrence, which, from their nature, produce 
expressions, each invariant under orthogonal transformation in the sense that 
after the transformation has been performed on the various quaternions of 
which it is a function, the result is the transform of the initial function (§ 8). 
Certain of these forms are also invariant under any linear homogeneous trans- 
formation. Similar forms were used by MAcMAHONff in the study of the 
ternary and quaternary invariants and covariants under orthogonal trans- 
formations. The study of differential parameters and their syzygies is made 

* Presented to the Society (Chicago), April 8, 1910, and April 28, 1911. 
jAmerican Journal of Mathematics, vol. 19 (1897), pp. 193-216. 
t Elements of Quaternions, 2d ed., vol. I, page 563. 
§ These Transactions, vol. 2 (1901), pp. 183-214. 
|| These Transactions, vol. 3 (1902), pp. 46-59. 
If Manual of Quaternions, page 264. 

"Wiener Sit zungsbericht e, vol. 115 (1906), Abt. lla, pp. 739-775. 
ft Proceedings of the London Mathematical Society, ser. 2, vol. 1 
(1904), pp. 210-229. 
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simpler by the use of these forms, and the symbolic invariants of Maschke * 
appear as a natural part of the development. Two differential operators are 
introduced which are the extensions of the quaternion V to space of four di- 
mensions and to curvilinear three-dimensional space. 

Up to the present time, the greater part of the development of the quaternion 
algebra has been in terms of scalars and vectors rather than in terms of quater- 
nions proper. This often introduces prolixity in the development. It is hoped 
that the present paper will present some fruitful extensions along the more gen- 
eral line. In the expansion of the forms appearing here it will be evident that 
they are of the nature of determinants and pfaffians. That other forms 
might be studied is of course self-evident, but these are the ones of most 
frequent occurrence in the functions of geometry and physics. However the 
developments below are on the formal side and valid irrespective of the inter- 
pretation of the quaternions. 

Part I. Fundamental Formula. 

1. Throughout the paper italic letters will be used to represent quaternions, 
which may occasionally reduce to scalars or vectors. The conjugate of a 
quaternion q will usually be denoted by q' , but the conjugate of an expression, 
as qrs, will be denoted by K • qrs. 

2. We notice that the quaternions uv' and v' u are important, particularly 
on account of the relations shown in the formula? 

uv' = S • uv' + YuSv — YvSu — YYuYv , 

which resolves uv' into three parts whose products have vanishing scalars; 
and if I and m are scalars, 

uv' ( lu + mv ) = u ( 2lSuv' + viTh ) — ZiT 2 u = ( lu + mv ) v'u , 

which shows that any quaternion of the form lu + mv is converted into another 
of the same form by left-hand multiplication by uv' or by right-hand multi- 
plication by v' u . Also uv' and v' u differ only in their axes, having the same 
tensors and angles. The quaternions uv' and vu' however differ only in hav- 
ing opposite angles. If r = lu + mv , we have uv'r — rv' u = . 

3. We define first the two functions 

I • uv = \ ( uv' + vu' ) , A2 • uv = 5 ( uv' — vu' ) . 

The first of the two functions is obviously a scalar, since it is the sum of a 
quaternion and its conjugate, and the second is a vector. Indeed, 

u ?/ 

I • uv = S ■ uv' = TH ■ S • - , A 2 • uv = Yuv' = T*v ■ Y ■ - . 

V V 



* These Transactions, vol. 1 (1900), pp. 197-204; vol. 7 (1906), pp. 69-80. 
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From the definitions we have at once 

I • uv = I • vu = I • u' v' = I • v' u' , A2 • uv = — A2 • vu . 

The latter expression is alternating and therefore vanishes if u is a scalar 
multiple of v . If we write p = XJWuYv, we find with no difficulty that 

A2 • u' v' = p • A2 uv • p _1 . 

Hence A2 • uv differs from A2 u' v' only in direction. By expanding the two 
sides of the equations, we find that 

U~ X Ao • uv • u = A2 • v' u' , V -1 • A2 • UV • V = — A2 • u' v' . 

The two forms used by Joly are, in this notation, as follows : 

( uv ) = vSu — uSv = \ ( A2 • v! v' — As • uv ) , 
[uv] = YYuYv = \ ( — A 2 • u' v' — A 2 • uv) . 

The single form given by Hamilton, [uv], is twice the corresponding one 
given by Joly, that is, is 2 YYuVv. 

If I • uv = we shall speak of u and v as orthogonal. The solution of the 
two equations which express that r is orthogonal to both u and v is 

r = x ( uv' — v' u ) + y ( uv' uv' — uv' v' u ) , 

where x and y are arbitrary scalars; for, this value of r satisfies the two equa- 
tions, and r can depend only upon two arbitrary parameters. 
4. Passing now to three quaternions, we define the functions 

Ai • uvw = \ ( uv' w + wv' u ) , A3 • uvw = \ (uv' w — wv' u ) . 

It is evident that Ai • uvw = Ai • wvu. Also 

Ai • UVW = 2 ( vw> u ~H wv ' u— vu' w — vw' u + uv' w + vu' w) 

= u\ • vw — vl • wu + wl • uv. 

Again from the definition it is obvious that 

A3 • uviv = — A3 • wvu. 

Since we have identically 2wl • vw = 21 • vw ■ u, we have by expansion the 
identity 

uv' w + uw' v = vw' u + wv' u , 

from which by transposition, and by reapplying the result, we arrive at 

uv' w — wv' u = vw' u — uw' v = wu' v — vu' w . 
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Therefore, returning to the form A3 , we have 

A 3 • uvw = A 3 • vwu = A 3 • wuv = — A3 • wvu = — A3 • uwv 

= — A3 • vuw . 

This function is therefore an alternating function of u , v , and w . It vanishes 
if they are linearly connected. Joly's form [uvw] is in this notation 
A3 • u' v' w' . His form (uvw) is SA3 • u' v' w' . We may expand in the form 

A 3 • uvw = I [uv' w — uw' v + vw' u — vu' w + wu' v — wv' u ] 
= J [uAi • v' w' + vKi • w' u' + wki ■ u' v' ] . 

This form shows that the function is of the nature of a determinant.* Again, 

I • MA3 • uvw =0=1- vh.z • uvw = I • wAz • uvw. 

For, if w,e expand I • MA3 • uvw we have 

\ [u' uv' w — u' ivv' u-\-w'vu'u — u' vw' u ] 

= j [u' w.I • vw — u' I • vw ■ u] = . 

A similar proof holds for v and w . Hence the solution of the three equations 
expressing the orthogonality of r to u , v , w , is r = £ A3 • uvw , where x is an 
arbitrary scalar. 

We see also from these equations and the expanded form of Ai • uvw that 
we may resolve uv' w into two orthogonal parts : 

uv' w = Ai • uvw + A 3 • uvw . 

If Ai • uvw = , we have uv' w + wv' u = , whence 

uv' tvw' v + wv' uw' v = = wv' uv' w + wv' wv' u . 

Adding the two sides of this double equality we have 

[wv' wv' + wv' ( w' v + v' w ) + vv' ww' ] u = . 

This gives easily wv' u\ • vw = 0. Similarly we find that uv' wl • uw = 0. 
Hence either we must have I • vw = = I • uw or else uv' w = = W tt . 
The latter cannot hold unless T 2 u , or T 2 » , or T 2 w vanishes. Then if the 
quaternions are real we must have, when Ai • uvw = 0, the three equations 
I • uv = = I • uw = I • vw , and the three quaternions form an orthogonal 
system. 

We find thatf 

TAi • uviv = V ( ImmI 2 vw + IotI 2 wu + Iwwl 2 uv — 21uvlvwlwu ) 

= TAi • vwu = • • • . 

* Cf . Joly: Transactions of the Royal Irish Academy, 1902, 32A, Part 
II, pp. 17-30. 

jThe period following I and A will be omitted when no confusion can result, just as it is 
after S and V in treatises on Quaternions. 
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Since T 2 uv' w = 1 when u , v , and w are unit quaternions, and 

T 2 uv' w = T 2 Ai • uvw + T 2 A3 • uvw , 

we might call TAi • uvw the cosine of the triple u, v, w and TA3 • uvw the 
sine. The latter for geometric reasons could also be called the Staudtian. 
Moreover, 

luu luv \uw 

T 2 A3 • uvw = I • A3 • uvw A3 • uvw = Ivu Ivv \vw 

Iwu Iwv Iww 

5. If we write out each quaternion in terms of 1 , i , j , k we find that 

1 i j k 



A3 • uvw = 



Uq U\ U<l Uz 
Vq Vi Vi V 3 
Wo Wi Wz w 3 



From this we have at once 

I • £A.3 • UVW = \t Ul » 2 «>3 I . 

Also from the properties of determinant multiplication we have the important 
formula 

lau lav law 



I • A 3 • abcAz • uvw = 
We may also easily verify that 

A3 • bcAz ■ uviv = 



Ibu Ibv Ibw 
leu lev lew 



— u — v — w 

Ibu Ibv Ibiv 

leu lev lew 



Hence we have at once 

I • aA 3 bcAz uvw = — IA 3 abcAz uvw . 
6. Some useful formulae are 

A 2 • A 2 uv A 2 wx = — A 2 uwlvx + A 2 uxlvio + A 2 vwlux — A 2 vxluw , 

Ai • uAz vwx = — A 2 vwlux — A 2 wxluv — A 2 xvluw , 

A 3 uvw • u' = mA 3 u' v' w' , 

A 3 ■ A 2 uvAi vwAi wu = T 2 A 3 uvw . 
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7. For four quaternions we define the three functions 

I • uvwx = IwAi vwx = luvlivx — luwlvx + luxlvw 

= \ ( uv' wx' + ux' wv' + vw' xu' + vu' xw' ) = \ ( S • uv' wx' + S • vw' xu' ) . 

A 4 • uvwx = — luAs vwx = I ( uv' wx' — ux' wv' + xiv' vu' — vw' xu' ) 

= \ ( S • uv' ivx' — S • vw' xu' ) . 

A2 • uvwx = \ (uv' wx' — xw' vu' ) = V • uv' wx' . 

Evidently in luvwx we may permute the four quaternions cyclically, or may 
reverse the order. A4 is an alternating function. In A 2 the order may be 
reversed by changing the sign. luvwx is a pfaffian, for it is the square root of a 
skew symmetric determinant of even order; indeed we have by actually squar- 
ing the expanded form, 



( luvwx ) 2 = 






luv 


luw 


— Ivu 





Ivw 


— Iwu 


— Iwv 





— Ixu 


— Ixv 


— Ixw 



lux 

Ivx 

Iwx 





8. If now we define a function A5 • uvwxy in the same manner as above it will 
be alternating and vanish, since any five quaternions are linearly connected. 
Indeed the vanishing of this expression gives us the important identity 

mA 4 vwxy — 0A4 uwxy + wA 4 uvxy — aA 4 uvwy + yA$uvwx = . 

Operating by I • z( ) , where z is any quaternion, we find the important identity 

zA 4 vwxy = — A 3 wxylvz + A 3 vxylwz — A 3 vivylxz + A3 vwxlyz . 

We define further 

Ai • uvwxy = 2 ( uv ' wx> V + V x> wv' u ) > 

A3 • uvwxy = 2 ( uv ' wx ' y ~ y x> wv ' u ) • 
Functions like these and the even numbered functions I, A2, A 4 may be 
defined for any number of quaternions, and are useful in reductions. 

As previously stated, all these forms are invariant or pseudo-invariant 
under the substitution of aqb for q, where Ta = Tb = 1. This is evident 
by mere substitution in the definitions. For example , 

A2 uv = a' A2 (aub) (avb) • a, A\uvw = cT 1 ■ K\(aub) (avb) (awb) • b" 1 , 
A3 • uvw = a' • A3 ( aub ) ( avb ) ( atvb ) • b' , 

while I • uv , I ■ uvwx, A 4 • uvwx are invariant. 

The functions A3 uvw and A 4 uvwx are invariant under any linear homogene- 
ous substitution, in the sense that such a substitution merely multiplies the 
function by the determinant of the substitution. 
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Part II. Differential Operators. 

9. Let us suppose that the coordinates of a point in space of four dimensions 
are w, x, y, z, corresponding to the real or scalar axis and the vector axes 
i , j , k . We then define the differentiating operator * D as follows : 

d . d d d 

D = - — h«7+jr + *v. 
dw dx dy dz 

In physical applications it would be desirable to define this operator directly 
in physical terms, and in geometry in geometrical terms, but for the present 
formal developments this definition seems simpler. It is obvious that if the 
quaternion q becomes q + dq , where 

dq = dw -\- idx + jdy + hdz , 
then the operator 

I • dqD = dw— + dx — + dy — + dz — 
dw dx dy dz 

gives the differential change in Q , any function of the position q , due to the 
differential change dq. If Q is a scalar function of position and its levels 
given by Q = c , then the change in Q which is greatest will be normal to these 
levels and will be given by DQ itself. Let 

Q= W + iX + jY+kZ. 
Then 

1 V V dw + dx + dy + dz • 

We might call this the divergence of the function Q . Its vanishing gives an 
equation well known in hydrodynamics, the equation of continuity. In this 
case we take Q = p ( 1 + <r ) where a is the velocity, p the density, w the time. 
If we combine D with itself we have 

*\2 r*2 n2 n2 r\2 n2 n2 J}2 

I ' DD= dw i+ dx~ 2 + ~d^' i+ dJ 2 ' I - DD ' = -Q—2- -Q-T- frf--^- 

The first is the extension of the Laplacian to four dimensions and leads to what 
may be called four-dimensional potential functions. W T e have an example of 
the vanishing of these differential forms in the differential equation of wave 
motion, which if we set w = at may be written 

I • DD • u = 0, or I • DD' • u = 0, 

according as the term containing the differential of the time is negative or 
positive. 

* This is the operator lor of Minkowski, Mathematische Annalen, vol. 68 
(1909-10), pp. 472-525. 
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The first partial derivatives of a scalar function Q as to the coordinates, 
are the coordinates of the quaternions 1 , i , j , k in DQ . The second partial 
derivatives would be the coefficients of the dyads i n DIDQ . In fact, the 
scalar invariants f of the linear quaternion differential operator DI • D ( ) , 
namely Mi , M 2 , M a , M 4 , lead not only to the Laplacian, but to other forms 
of operators, all of which are invariant under orthogonal transformations 
of the axes 1 , i , j , k . 

We have 

T^n-f^K a _ZiiZ M\i -( bW Hi dZ ?A\ 

A °- DQ ~ l \ dx ~ dw + dz ~ dy ) + J V By ~ dw + dx~ dz ) 

(aw dZ dX dY\ 
+ \ dz dw + dy dx ) ' 
a ivrv ( BW^dX^dY dZ\,.( dW^dY^dZ dX\ 
A ^ Q=l {—dx~ + ^+-dz-'dy) + A-^i + ^ + ^~^) 

( dW dZ dX dY\ 
+ \ dz + dw + dy dx ) " 

We may say that if A2 DQ = , Q is leftward irrotational, if A2 D' Q' = , 
then Q is rightward irrotational. The expressions A2 we may speak of as 
left curl and right curl. If both vanish we may speak of Q as irrotational. 
The rotations here are of course planar and not axial. As an example of the 
use of these forms, let Q = <p — A be a quaternion potential function, <p being 
a scalar potential, and A a vector potential. Then 

A 2 DQ = v^ + dA I dw + VvA = - x + y , 

A 2 D' Q' = - V<p - dA I dw + VvA = x + y . 
Thence 

2x= A 2 D' Q' - A 2 • BQ , 2y = A 2 D' Q' + A 2 T)Q . 

From the identity A 3 DQD = we have, taking scalars, S • VDVQVD = , 
or 

Svy = . 
Taking vectors we have 

V (A 2 DQ • D + DA 2 D' Q') = 0, VD' (- x + y) + VD (a; + y) = 0, 
that is, 

dw 

Both equations are thus included in A 3 DQD = 0, an identity. 

Setting 2DA 2 • Q'D' = A 2 DQ • D - DA 2 D' Q' = 2q, we have, by taking 
scalars, 

SVx = — ■ Sq. 

t These invariants are coefficients of the characteristic equation; cf. Joly, Manual of Qua- 
ternions. 
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Taking vectors we have 

These two equations arise thus from a single equation. 

Now if we set in these f ormulse w = V — 1 • ct , and in the first place set 
x = — V — l-E,y = H,Sq=V / — 1 • p , we have the common equations 
of the electro-magnetic field.* If we set in the next place x = — V — 1 • E , 
y = B,ovx=H,y= — V — 1 • D , respectively, we arrive at the vectors 
called by Minkowski f electric force of rest, and magnetic force of rest, 

E = £(A 2 DQ.r-rA 2 D'Q'), S = A 3 DQr, 
provided r is the velocity given by 

V - 1 • r = (1- V -1 • q)/ V (1- q 2 ). 
If we start with such equations as the familiar ones 

Vvtf- — =p<t, Sv£=-p, YvE + ~ = 0, Svtf=0, 
we may reduce them to the present forms. We have first 

TI dE 77 , 9H 

VH- Tt = p «, -Ev + - dt = P . 

Adding these, we have the single equivalent equation DH — ED = q . If 
// = YvA , and E = W + dA / dt this reduces with little trouble to 

D'ID'Q- IDD' • Q' = 2q. 

Formula for the extensions of Stokes' theorem and Green's theorem will 
be found in the reference to Joly below.J 

When we operate on Tq~ 2 or ( qq' ) _1 the operator IaD will give the hyper- 
spherical harmonics, and the results are susceptible of interpretations similar 
to those for three dimensions. The introduction of doublets follows analogous 
lines. These applications can only be mentioned. § 

* Cf . Lewis, Proceedings of the American Academy of Arts and 
Sciences, vol. 46 (1910), pp. 165-181; M. Abraham, Rendiconti del Circolo 
Matematico di Palermo, vol. 30 (1900), pp. 33-46. A. Sommerfeld, Annalen der 
Physik und Chemie, ser. 4, vol. 32 (1910), pp. 749-776; vol. 33, pp. 649-689. 

t Minkowski, Gottinger Nachrichten, 1908, pp. 53, 111. The imaginary is not 
really necessary. 

JProceedings of the Royal Irish Academy, ser. 3, vol. 8 (1902), pp. 6-20. 

§ Cf . Joly, Appendix to Hamilton's Elements of Quaternions, 2d ed., Appendix 13. Tait, 
Scientific Papers, vol. 1, 37-42, 134-135, 136-150, 176-193, 234-236, 280-281; vol. 2, p. 312; 
Quarterly Journal of Mathematics, vol. 6 (1864), pp. 279-301. McAtjlay, 
Messenger of Mathematics, vol. 14 (1884), pp. 26-37; Proceedings of the 
Royal Society of Edinburgh, vol. 18 (1890), pp. 98-123. 
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10. We pass now to the development of certain differential formulae for 
three-dimensional curvilinear space, or what is the same thing, quaternions 
depending upon three parameters. This is equivalent to considering the 
three-dimensional space in a flat four-dimensional space, so that the problem 
becomes that of the differential geometry of curved three-spreads in four- 
dimensional space.* Let us suppose, then, that 

q = /Oi, «2, u s ). 

This equation will limit q to a certain three spread. We shall write 

tf tf tf i , , , , 

qi = —— , q 2 = -r— , q 3 = -r— , dq = qi du x + q 2 du 2 + q 3 du 3 . 
dtli du 2 du 3 

The quaternions q\ , q 2 , q 3 are independent, as the parameters are essential. 
Let the quaternion A 3 qi q 2 q 3 = n. Then n is the normal at q to the three- 
flat tangent to the three-spread, since Inqi = , Irag 2 = , Inq 3 = . The 
differential dq lies in the tangent three-flat. 

11. We shall indicate the quaternions 

— A 3 nq-i q z — A 3 nq$ qi — A 3 nqi q 2 

Inn ' Inn ' Inn 

respectively by <ji , q 2 , and q 3 . These give at once 

Iqiqi=l, Iqtqj = (i #= j) , Inq t = . 

Any quaternion in the three-dimensional space may be expanded in the forms 

r = q\Iq\r+ q 2 Iq 2 r + q 3 Iq 3 r = qi Iqi r +q 2 Iq 2 r + q 3 Iq 3 r . 

We now define the operator 

We have by substitution of values the operator 

IdqA. = dui \- du 2 \- du 3 —— . 

dui du 2 du 3 

That is to say, the operator IaA gives the rate of variation of any operand due 
to a variation in the three-spread of q in the direction a, where a is a unit 
quaternion. It is evident again that if Q is a point function defined for the 
three-spread, a scalar let us say at first, then it will have certain level surfaces. 
The quaternion AQ gives us then the rate of maximum change, which will be 
the normal to a level surface. The operator A is thus the extension to curved 
* Thus, if q is to terminate on a hypersphere, we might write 

q = u x i + u 2 j + « 3 k + v ( a 2 — u\ — u\ — u\ ) . 
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space of the operator V for flat three-dimensional space. We proceed to show 
that it is invariant under transformations of the parameters. To do this let 

Ml = <p! ( Vi , Vi , V 3 ) , M 2 = <P2 ( n , V 2 , V S ) , Us = (p 3 ( V! , Vi , V 3 ) , 

where the functions tp are differentiable and without singularities at the points 
considered. Then 

A ^A d(v 1; V 2 ) ( dq \ 

A29298 = SA291 " 92 ' •*(«„«,)• '" \*'-a^-)- 
Thence building up A v we have A„ = A„ . Also we have 

Vn„ = UA391 q 2 q 3 = UA 3 (qi) v {qi)v (qa)v = Un„ . 

The normal is thus an invariant under transformation of the parameters. 
It follows naturally that such linear operators as I ( ) A • Un , I ( ) A • Ax 
and powers and combinations of these are also invariant. As an example, we 
see that the function IaA • Q gives the rate of variation of Q at the point 
in the direction o . The axes of this function are the directions of extremal 
variation and the roots are the rates of such change, at least if the function 
is self-transverse. If we lay off in each direction the value of the function 
( when Q is a scalar) the points so determined will have a hyper-quadric for 
their locus. The method of treatment differs little from the corresponding 
case for three dimensions. 

Again consider the function N = I ( ) A • Vn . This is evidently the rate 
of change of the unit normal at the point. From the character of A it 
is self-transverse. Hence the axes are orthogonal and are the directions of 
the curvatures (reciprocals of the radii) at the point. The roots are real and 
are the curvatures themselves for real quaternions gi , q 2 , q 3 . One root 
is zero, the corresponding axis being the normal. The usual scalar invariants 
of N are therefore the mean curvature, the second mean curvature, and the 
total or Kronecker-Gaussian curvature. The mean curvature (as in the 
similar three-dimensional case) is IAUm. The second mean curvature and 
the total curvature are respectively 

I • A 3 nA' A" A 3 nr 1 LW LV , IA 3 A' A" A'" A 3 LV Um" Vn'" , 

where the accents are dropped after the expansion. 

12. Since Am, Av, Aw, ••• are invariant, also TAm, TA«, ••• are invari- 
ant. Indeed it is evident that the differential parameters of the first order 
for the functions u , v , • • • are, in the usual notation for differential para- 
meters, 

Ai u = IAmAm, Ai (m, v) = lAuAv. 

In addition to these we have the following extensions of the Darboux ( u , v) 
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for two dimensions to three dimensions : 

A 2 AuAv , Ai AuAvAw , A 3 AuAvAw, 

giving the respective parameters 

IA 2 AmA»A 2 AuAv , I • Ai AmAdAwAi AuAvAw , IA 3 AmA»AmA 3 AuAvAw , 

IAuAvAwAx , IA 2 AuAvAivAxAz AwAsAwAa; . 

The form A 4 AuAvAwAx vanishes identically. 

The formulae of part I give various syzygies connecting these, for example 

I • AuAvAwAx =Ai(w, v) Ai (w , x) — Ai (m, w) Ai (v , x) 

+ Ai (m, x) A x (v, w) . 

Since Am is the normal in the three-spread for the levels of the scalar function 
u , the parameter of first order of u , Ai m , is the square of the length of this 
normal. The differential parameter of first order of u and v is the projection 
of the normal to the level of either, on the other, multiplied by the length of 
the latter. If it vanishes the levels are orthogonal. The interpretation of 
the other forms follows the exactly analogous forms for V in a flat three- 
dimensional space. 

13. Any parameter of the first order is itself subject to operation and we 
arrive thus at compound forms such as lAuA ( Ai u) . We are led thus to 
consider operators such as $ = I ( ) A • Am . The first scalar invariant of 
this linear quaternion operator is lAAu , which is the well-known differential 
parameter A2 u . It vanishes if u is a solution of the generalized Laplacian 
equation. Indeed if the space becomes flat and the position is determined by 
p instead of q , the second differential parameter becomes — V 2 . If the lines 
u\ , it* , U3 are isothermal, we have a simplification which is not difficult to 
follow. The second invariant of $ is 

I A 3 nA' A" A 3 rT 1 ( Am)' ( Am)" . 

This can be expanded in the usual manner for the coefficients of the charac- 
teristic equation of a matrix or linear quaternion operator, in the form 

\ (IAAm) 2 - IA' A" mIA' A" m . 

It is the well-known parameter A22 m . Finally the third invariant is 

- £ IA 3 A' A" A'" A 3 (Am)' (Am)" (Am)'". 

This might be called the parameter A222 m . The fourth invariant vanishes 
identically. 

The function $ is the rate of change in any direction of the vector rate 
of change of the function which is itself the normal to a level of the function. 
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14. If we consider only the direction of the normal of the level, we have the 
function E = I ( ) A • UAw , which gives the rate of deviation of the normal 
direction. There are now only two invariants that do not vanish identi- 
cally, and these determine what may be called geodetic curvatures. The 
first invariant is iAUAw , the geodetic mean curvature. The second is the 
expression 

IA 3 nA' A" A 3 n' 1 UAm' UAw" . 

These are the extensions of the usual geodetic curvature for two dimensions. 

15. We notice next the fundamental forms or co variants of different orders. 
The first is 

Idqdq = 2lgi q\ du\ dui + 22l</i q» du\ du 2 . 

The coefficients Ig; qj correspond to the E , F , G of the geometry of sur- 
faces. They are expressible as parameters, for example, 

Igi q 2 = IA 3 nAu 2 Aw 3 A 3 nAw 3 Ami • 

Other covariant forms are 

Idqdr , A 2 dqdr , Ai dqdrds , A 3 dqdrds , Idqdrdsdt , A 2 dqdrdsdt , 

A 2 ndq , Ai ndqdr , A 3 ndqdr , A 2 ndqdrds , A4 ndqdrds = . 

Denoting differentiation as to Ui by the subscript i we have, from 
IUnUra = 1 , 

l\Jn(\Jn)i= 0, 

IUra (Ura)i = 0, IUra (Un) y = - I (Un)< (Ura)y (i,j = 1, 2, 3). 

Since Inqt = 0, therefore 

Inqa = — Iraq, = — Irijqi. 

In this we might also have written Vn instead of n . We arrive thus at the 
extensions of the D , D' , D" , of two dimensions. The Christoffel symbols are 

THI T [HI , T _ 

\ I \ = miki, \ i\ = 1 ?" AM i = iQkiqi- 

( ikrs ) = IA 3 rT 1 UniUn/cAznqrqs . 

Differentiating these and remembering that also A is subject to differentiation, 
we arrive at the symbols of higher order. This is the covariantive differenti- 
ation of Maschke.* 

The second fundamental form is 

IdqdUn = IdqNdq = Slgi ( Un ) 1 du\ du\ + 22lg, ( Un )j dui duj . 



* Maschke, loc. cit. Also these Transactions, vol. 4 (1903), pp. 457. 
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This form vanishes for the asymptotic directions, or the principal tangent lines 
of the three-spread. It is worth noting that the differential equation of the 
lines of curvature is A2 dqNdq = . This corresponds to the similar form 
for surfaces \dp<pdp = . 

16. The symbolic invariants of Maschke are easily expressed in these forms 
by noticing the following equivalences * for his symbols 

fi = Uqi, fifk = Ifgiir<Zfc = iqiQk, fn = lf?«. 

(/) = A 4 Unfits, 
whence 

(/) 2 = 3 ! ItMJw = 3 ! , (uf ) = A 4 nAwf 1 f 2 . 

(fu ) = A 4 VnAu' Au" f 1 , 
whence 

(uff = 2lAuAu = 2AiM. 

( (uf) , /) = A 4 UwA (A 4 UnAwfi f 2 ) fi r 2 = 2IUwUnlAAM = 2A 2 m . 

Not to elaborate too far, we need only to remember that each ( ) is an A 4 , that 
u, v , w , x , etc., are functions of q such as we have been using, that for every 
/ or like symbol we substitute Ifg, for subscripts we differentiate as to the 
corresponding u. All the Maschke symbols for four dimensions become at 
once interpretable in this system. The vectors of Ingold f are the quaternions 
appearing in these forms, his formula? being analogues of Maschke's, not 
equivalents. The ordinary formulae hold of course for these quaternions, but 
the Maschke forms are all scalars and commutative, and are the scalar forms 
with the f-pairs introduced. 

* The bilinear function Q ( £ , f ) is an abbreviation for Q ( <?i , </i ) + Q ( g 2 , g 2 ) + Q ( g 3 , 93 ) . 
We may always write for a bilinear function Q ( f , f ) the form Q ( A, q ) or equally 
Q (q, A ) . See McAtjlay, Utility of Quaternions in Physics ; Shaw, Synopsis of Linear 
Associative Algebra. 

t These Transactions, vol. 11 (1910), pp. 449-474. 



